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Abstract. The differential cross sections cro — o-t+£ctl,o-lt, and <jtt of it electroproduction from the proton were 
measured from threshold up to an additional center of mass energy of 40 MeV, at a value of the photon four-momentum 
transfer of Q 2 = 0.05 GeV 2 /c 2 and a center of mass angle of 9 = 90°. By an additional out-of-plane measurement 
with polarized electrons a LT i was determined. This showed for the first time the cusp effect above the ix + threshold 
in the imaginary part of the s-wave. The predictions of Heavy Baryon Chiral Perturbation Theory are in disagreement 
with these data. On the other hand, the data are somewhat better predicted by the MAID phenomenological model and 
are in good agreement with the dynamical model DMT. 
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> ! 1 Introduction 

5— I Over the last two decades, neutral pion photo- and electropro- 
. 5^ i duction at threshold have been established as a testing ground 
for Chiral Perturbation Theory. First experimental attempts at 
Saclay [ 1 1 and Mainz Q to test the low energy theorem for the 
s-wave multipoles at threshold [30 initiated the development 
of an intense investigation of the threshold region with close 
cooperation between experiment and theory. 

Chiral Perturbation Theory (ChPT) is the effective field the- 
ory (low energy approximation) of QCD [5,6,7,8]. The first re- 
liable ChPT calculations of pion photoproduction in the direct 
threshold region showed that the old low energy theorem has to 
be modified, and that the s-wave prediction is a slow converg- 
ing series in this calculation [9]. For two p-wave combinations, 
new low energy theorems were derived [ 10 1. 
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Fig. 1. The rescattering process responsible for the cusp in the 7r° 
s-wave amplitude at threshold. 



The energy dependence of the s-wave multipoles in the 
7P — > ir°p reaction has a smooth direct part and a rapidly vary- 
ing (cusp) part due to a ir + rescattering contribution shown in 
Fig. [U[ 1 l ]fl2lfT3ll . Since the ir + production amplitude is domi- 
nated by the additional coupling of the photon to the charge of 
the pion, this amplitude is much larger than the production am- 
plitude of neutral pions, leading to a sizable unitarity cusp in 
the 7T° s-wave amplitude which appears in the imaginary (real) 
part above (below) mr + threshold. 

The predicted cusp effect in the real part of the transverse 
s-wave amplitude Eq + , was established in photoproduction ex- 
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Fig. 2. Definition of the pion angles 9 and (f> in the laboratory frame. 



periments at SAL OH and MAMI lll5lll6lfT7l . While these 
experiments were only sensitive to the real part of this am- 
plitude, one aim of this work is to access the imaginary part 
of the s-wave amplitudes. In electroproduction, this can be 
achieved by an out-of -plane measurement with a longitudinally 
polarized electron beam to access the interference cross sec- 
tion ctlts which is proportional to the imaginary part of the 
transverse and longitudinal s-wave amplitudes Eq + and Lo + , 
respectively. 

The second aim of the present work is to improve the ex- 
perimental determination of the p-wave amplitudes. The exist- 
ing electroproduction experiments with large solid angle lfT8l 
\l9l covered the energy range only up to the first 4 MeV above 
threshold. By an extension of the experiment to an energy re- 
gion of 40 MeV above threshold, where the p-waves clearly 
dominate, a more accurate separation of the cross sections is 
possible. 



systematic errors: 

LT ' { ' G++a- a T (9) + ea L (9)-ea TT (e) (J 

where a + and a~ are the differential cross sections for <\> — 
90° with beam polarization parallel and antiparallel to the beam 
direction, respectively. 

In the threshold region, all multipole amplitudes with an- 
gular momentum I > 1 are negligible. The cross sections can 
be further decomposed in s- and p-wave multipoles as 

a T {9) = p/k 7 (A + B cos 9 + C cos 2 9) , 
cj L {9) = p/k 7 (A' + B' cos 9 + C cos 2 9) , 
cr TT ((9) =p/fc 7 (Fsin 2 9) , 
a LT (9) = p/fc 7 (D sin 9 + E sin 9 cos 6) , 
a LT , (9) = p/k 7 (G sin 9 + H sin 9 cos 9) , (4) 

where p is the pion center of mass momentum and fc 7 = 
(W 2 — m 2 + q 2 )/(2W) is the equivalent real photon en- 
ergy. The angular coefficients contain two s-wave and five p- 
wave multipole combinations [20 21 1. By measuring the helic- 
ity asymmetry with an out of plane measurement, we obtain 
a lt' which is proportional to the imaginary part of (Eq + P^ + 

In this energy region the p-wave multipoles are approxi- 
mately real so we are directly sensitive to the imaginary part of 
the s-wave multipoles. 



2 Notation and formalism 



3 Experiment 



In the one photon exchange approximation the electroproduc- 
tion cross section of pions can be written as (see e.g. |20|) 



da {6, 



= r (a T (9) +ea L {9)+e a TT {9) cos 20 
&lt(9) cos</> 



h^2e(l - e) a LT > (9) sin 



(1) 



with the virtual photon flux r, the transverse photon polariza- 
tion e, and the longitudinal beam polarization h. In addition, the 
cross sections also depend on the center-of-mass energy W (or 
AW = W — Wthreshoid) and on the photon four-momentum 



transfer q 



Q 



q , with the photon laboratory 



energy and momentum lj and q. 

The pion emission angles 9 and <f> are defined relative to the 
direction of the momentum transfer q. Figure|2]shows the def- 
inition of the angles in the laboratory frame. In the following, 
these angles are used in the photon-proton cm. system. 

In this experiment, the photon polarization e was not varied, 
so the differential cross section 



a (9) = a T {9) + ea L {9) 



(2) 



could not be separated further. 

Since ctlt' is small compared to the unpolarized cross sec- 
tion, it is useful to define the asymmetry Alt' {&) to minimize 



The experiment took place at the three-spectrometer setup of 
the Al Collaboration at MAMI (see E2J for a detailed de- 
scription). It was performed and analyzed based on the same 
techniques described in refs. [18, 1 9 1 . From Eq. [TJ one can 
see that for a separation of the cross sections, a measurement 
at <fi — 0°, 90°, 180° is the best choice. A measurement at 
4> = 90° requires out-of-plane detection. 

Spectrometer B of the Al setup can be moved out of plane 
up to a cartesian angle of 10° in the laboratory frame [23|. By 
choosing spectrometer B for the electron detection at a small 
forward angle, we were able to cover <fi — 90° up to a center-of- 
mass momentum of lOOMeV/c within three kinematical set- 
tings. In addition, four in-plane settings were chosen to sepa- 
rate the cross sections for = 0° and <f> — 180°. 

Table [TJ summarizes the kinematical settings including the 
approximative acceptance range in the center-of-mass momen- 
tum. For all settings, the four-momentum transfer of the photon 
was fixed to the value of Q 2 = 0.05 GeV 2 /c 2 used in former 
experiments 1191 . At the beam energy of Eq = 854.49 MeV, 
adjusted to optimize the spin precession angle between polar- 
ized beam source and spectrometer setup, the photon polariza- 
tion was e = 0.933. 

The polarized electron source |24| of MAMI was run at a 
polarization of 75% at an average current of 5 fiA. The beam 
polarization was determined by an elastic H(e, e'p) measure- 
ment with a focal plane polarimeter in spectrometer A to an ac- 
curacy of ±5% absolute. To avoid systematic errors, the beam 
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Table 1. Kinematical settings for the central trajectory of the spec- 
trometers. The symbols 6 e and <j> e mark the in-plane and out-of- 
plane cartesian angles of spectrometer B. Common for all settings are 



Q 2 = 0.05 GeV 2 /c 2 , 9 = 90°, e = 0.933, E b . 



854.49 MeV. 
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helicity was switched stochastically with a mean frequency of 
1 Hz. 

A subcooled liquid Hydrogen target with a oblong shape 
and a length of 5 cm was used. The beam was moved across the 
target in a controlled way by a fast raster magnet to avoid boil- 
ing of the Hydrogen. The beam current of 5 fiA corresponds to 
a luminosity of L = 6.7 MHz//.tbarn. A total effective beam 
time of 400 h was accumulated. 

Spectrometer A with a solid angle acceptance of 2 1 msr 
was used for the detection of the recoil proton, and spectrome- 
ter B with an acceptance of 5.6 msr was used for the detection 
of the scattered electron. The focal planes of both spectrome- 
ters were equipped with four layers of vertical drift chambers 
for position and angular resolution and two layers of scintilla- 
tors for trigger and timing resolution. In addition, spectrome- 
ter B was equipped with a gas Cerenkov detector for the sup- 
pression of charged pions. An overall momentum resolution of 
Sp/p < 10~ 4 and an angular resolution of better than 3 mrad 
was achieved. After correction for the flight path in the spec- 
trometers, a coincidence time resolution of 2 ns FWHM was 
determined. 
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Fig. 3. Reaction identification by coincidence time (a) and background 
subtracted missing mass (b). The events in the dark shaded area were 
defined as true coincidences, while the events in the light shaded area 
were used to subtract the background of random coincidences. 



4 Analysis and error estimation 

In a first step, the electron and proton were identified by the 
time of flight. The coincidence time was corrected by the flight 
path length inside the spectrometers of about 10 m. Figure[3](a) 
shows the corrected coincidence time distribution. The central 
peak of this distribution was used as the yield of true coinci- 
dences, while the events in the side bands were used to estimate 
the background of random coincidences. 

The reaction H(e, e'p)ir° was identified by a cut on the 
missing mass given by the four-momentum balance to^iss = 
(ein +Pin - e ou t -Pout) 2 - Figure[3](b) shows the missing-mass 
distribution after subtraction of the background of random co- 
incidences. The pion peak shows a missing-mass resolution of 
2.2 MeV/c 2 and a radiative tail to higher masses. The observed 
increase in radiative events at high masses was caused by the 
increase in acceptance and phase space. 

The integral of the accepted phase space was calculated 
with the Monte Carlo method. In this calculation, all resolution 



and efficiency effects and radiative effects according to the for- 
mulas of ref. Il27ll were included. The differential cross section 
a(9, <fi) was extracted for several bins in W, 9, and cf>. Since 
the acceptance in Q 2 and e was narrow, no binning in these 
variables was required. 

To separate the cross sections a x 2 -fit of the (^-dependence 
as given by Eq. [T]was performed for each bin in the center of 
mass energy W and with a cut on the pion angle of 75° < 8 < 
105°. 

By the cut in 9 a large fraction of events was lost. In ad- 
dition, the acceptance in 9 was different for each kinematical 
setting. To overcome these two problems, a second, model de- 
pendent method was used to separate the cross sections. In this 
method, the phenomenological model MAID ll25ll was used as 
parameterization of the cross section. For each event, the dif- 
ferential cross section was projected to the nominal kinematics 
at = 90° by 



cr(9O°,0) = cr(0,( 



MAID(90°, 
' MAID(6,q 



(5) 
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(a) Model dependent (b) Model independent 




10 20 30 40 10 20 30 40 



AW [MeV] AW [MeV] 

Fig. 4. The separated cross sections ctq, olt, and ott at 6 = 90°. Column (a) shows the result with model dependent analysis using the 
full statistics, while column (b) shows the model independent analysis with narrow kinematical cuts. The solid line shows the calculation in 
HBChPT (211, the dashed line shows the phenomenological model MAID (2003 fit) (25], and the dash-dotted line shows the DMT modelt26l. 
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By this method, the statistical error could be reduced for the 
price of an additional model error. 

The asymmetry Alt' was extracted by the helicity asym- 
metry defined in Eq. [3] including an additional background 
term in the denominator. The per-event-helicity and the ex- 
tended <f> acceptance were taken into account by a weight factor 
h sin (f> for each event. 

The systematic errors in detection efficiency, luminosity 
calculation, and reaction identification could be neglected com- 
pared to the statistical error. The only notable systematic error 
was assigned to the first three energy bins. Directly at thresh- 
old, the extracted cross section is very sensitive to the abso- 
lute momentum calibration of the electron arm, which was cal- 
ibrated to 300keV/c by a measurement of elastic electron- 
proton scattering. For the first three bins, this results in sys- 
tematic errors of 43%, 21%, and 2% respectively. For AW > 
3 MeV, the corresponding error is below 1%. For all other sys- 
tematic effects an overall systematic error of 3% was deter- 
mined 1221 . 

The model error of the model dependent separation of the 
cross sections was estimated by varying the multipoles of the 
model up to I = 1 separately by ±10%. This results in a model 
error in the extracted cross section of less than 5%. For the 
helicity asymmetry Alt' the dominant systematic error was 
caused by the measurement of the beam helicity h, which was 
known to 5%. This error is negligible compared to the statisti- 
cal error. 



5 Results 

Figure. |4] shows the separated unpolarized cross sections (To, 
<tlt, and ott- Column (a) shows the model dependent anal- 
ysis, column (b) shows the model independent analysis. Both 
results are clearly consistent within the error bars. The statis- 
tical advantage of the better utilization of the full data sample 
is compensated by the large systematic model error. Of course, 
the size of this error is disputable: the choice of a variation of 
10% in the leading multipoles is somewhat arbitrary and the 
authors are convinced, that this choice is conservative. In the 
future, better knowledge of all multipoles can clearly improve 
this analysis. 

Since the model dependent analysis shows that the model 
independent analysis is not affected by the acceptance match- 
ing, we will discuss only the model independent analysis. The 
data points of this analysis are compiled in Table|2] 

The cross sections are compared with three theoretical cal- 
culations. The solid line in Fig.|4]shows a calculation in Heavy 
Baryon Perturbation Theory by V. Bernard et al. ETI . As al- 
ready noticed in ref. lfl9l . this calculation overestimates the 
cross section at Q 2 = 0.05 GeV 2 /c 2 . The reason for this might 
be that some of the low energy parameters of this calculation 
were fitted to data at a photon virtuality of Q 2 = 0.1 GeV 2 /c 2 . 
This value may be somewhat too large for the validity of chiral 
perturbation theory. Therefore, a refit of the low energy param- 
eters to low Q 2 data might resolve some of this discrepancy. 
This would still leave the issue of the Q 2 dependence [ 19] as an 
unsolved problem. The more likely solution might be to carry 
out the ChPT calculation to higher order. The presently pub- 
lished calculation was carried out to order 0(p 3 ) for the p- 
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Fig. 5. The beam helicity asymmetry Alt' vs - center of mass energy. 
The solid line shows the calculation in HBChPT |21 1, the dashed line 
shows the phenomenological model MAID [ 25 j, the dash-dotted line 
shows the DMT model| 26|. 



waves with some extra counter terms of higher order for the 
s-waves. 

The two phenomenological models MAID [25 1 and DMT 
||26l both agree well with the unpolarized cross sections cfq and 
<jlt- The small interference cross section ott, however, is de- 
scribed only by the DMT model. 

Tableland Fig.[5]give the results for the helicity asymme- 
try Alt' ■ These data are the first experimental confirmation of 
the unitary cusp in an observable which is proportional to the 
imaginary part of the s-wave amplitudes. Unfortunately, due 
to the small size of this asymmetry the energy bins for this 
variable had to be made larger than for the unpolarized observ- 
ables and a detailed determination of the shape is precluded. 
The data point below nir + threshold is consistent with the ex- 
pected value of zero. Alt' is clearly underestimated by the 
HBChPT calculation ||2T1 . this is another hint that a higher or- 
der calculation is needed. The MAID model [25 1 overestimates 
Alt' by nearly a factor of 2. This clearly shows the necessity 
to include polarization observables into phenomenological fits 
to fix the small multipoles. The DMT model ll26l is in surpris- 
ing good agreement with the data. We want to stress, that this is 
a prediction without knowledge of the data. In this model, the 
threshold behavior is generated dynamically, taking the rescat- 
tering effects into account. 

In all of the three models | | 3> |Ps|, so that the primary 
sensitivity of A lt 1 is to the imaginary part of the longitudinal 
multipole Lq + . 



6 Summary 

In this experiment, we separated the cross sections oq, olt, 
and utt of neutral pion electroproduction at 9 = 90° at a 
photon virtuality of Q 2 = 0.05GeV 2 /c 2 from threshold up 
to 40 MeV above threshold. The cross sections were compared 
with predictions of Heavy Baryon Chiral Perturbation Theory 
||2T1 and the phenomenological models MAID [ 25 1 and DMT 
|. The discrepancy between HBChPT and the data suggests 



6 



M. Weis et al.: Separated cross sections in 7r° electroproduction at threshold at Q 2 — 0.05 GeV 2 /c 2 



the need for higher order calculations. The phenomenological 
models MAID and DMT show reasonable agreement with the 
separated cross sections. 

Simultaneous, the helicity asymmetry Alt 1 (90°) was mea- 
sured. This observable is sensitive to the imaginary part of the 
s-wave amplitudes, which was accessed for the first time by 
a direct measurement. This asymmetry can be described with 
impressive agreement by the DMT model. 
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Table 2. The separated cross sections ao, olt, and ott- Only the 
model-independent results with cut 75° < 6 < 105° are given in this 
table. 

AW <7 (90°) ctlt(90°) (7tt(90°) 

(MeV) (nb/sr) (nb/sr) (nb/sr) 



Table 3. The helicity asymmetry A LT > at 9 = 90°. 
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± 


29.4 


21.5 


806.0 


± 


35.5 


-245.5 


± 


21.1 


7.4 


± 


47.2 


22.5 


928.0 


± 


42.1 


-272.0 


± 


25.4 


-1.0 


± 


57.0 


23.5 


971.1 


± 


59.2 


-237.6 


± 


35.4 


-48.1 


± 


82.4 


24.5 


1018.3 


± 


46.2 


-256.4 


± 


27.6 


6.1 


± 


64.3 


25.5 


1113.4 


± 


79.0 


-250.8 


± 


48.9 


-70.6 


± 


113.6 


26.5 


1190.2 


± 


61.8 


-318.4 


± 


38.2 


24.0 


± 


88.9 


27.5 


1264.3 


± 


69.3 


-329.4 


± 


42.9 


-3.6 


± 


99.6 


28.5 


1329.0 


± 


79.2 


-331.9 


± 


47.2 


-16.9 


± 


113.4 


29.5 


1376.6 


± 


66.4 


-322.7 


± 


39.6 


-62.8 


± 


95.1 


30.5 


1427.1 


± 


75.0 


-318.1 


± 


44.7 


-135.2 


± 


107.4 


31.5 


1580.4 


± 


54.0 


-370.3 


± 


32.6 


-49.9 


± 


73.9 


32.5 


1722.7 


± 


54.7 


-422.6 


± 


31.4 


-121.1 


± 


73.8 


33.5 


1879.8 


± 


71.2 


-409.8 


± 


41.0 


-114.3 


± 


96.2 


34.5 


2003.2 


± 


70.2 


-470.2 


± 


40.3 


-167.6 


± 


94.7 


35.5 


2027.6 


± 


95.3 


-466.7 


± 


54.8 


-50.5 


± 


128.6 


36.5 


2169.0 


± 


110.7 


-464.3 


± 


63.7 


-209.0 


± 


149.4 


37.5 


2146.9 


± 


115.9 


-464.2 


± 


66.6 


-167.4 


± 


156.3 


38.5 


2223.2 


± 


127.1 


-477.1 


± 


72.1 


-51.4 


± 


168.3 


39.5 


2492.2 


± 


140.1 


-411.6 


± 


71.7 


-396.6 


± 


182.6 



